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Abstract 



The full set of equations governing the structure and the evolu- 
tion of self-gravitating spherically symmetric dissipative fluids with 
anisotropic stresses, is written down in terms of five scalar quantities 
obtained from the orthogonal splitting of the Riemann tensor, in the 
context of general relativity. It is shown that these scalars are di- 
rectly related to fundamental properties of the fluid distribution, such 
as: energy density, energy density inhomogeneity, local anisotropy 
of pressure, dissipative flux and the active gravitational mass. It is 
also shown that in the static case, all possible solutions to Einstein 
equations may be expressed explicitly through these scalars. Some 
solutions are exhibited to illustrate this point. 
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1 INTRODUCTION 



The gravitational collapse of massive stars, and its resulting remanent (neu- 
tron star or black hole) represent one of the few observable scenarios where 
general relativity is expected to play a relevant role. Therefore, a detailed 
description of gravitational collapse and the modelling of the structure of 
compact objects under a variety of conditions, remain among the most inter- 
esting problems that general relativity has to deal with. This fact explains 
the great attraction that these problems exert on the community of the rela- 
tivists. Starting with the seminal papers by Oppenheimer and Snyder [1] and 
Oppeheimer and Volkoff [2], a long list of works has been presented trying 
to provide models of evolving gravitating spheres and compact objects (just 
as a sample, and without the pretension of being exhaustive, see [3]-[26] and 
references therein). 

Motivated by the above arguments, we present in this work a study on 
self-gravitating relativistic fluids in terms of a set of scalars obtained from 
the orthogonal splitting of the Riemann tensor. As we shall see these scalars 
have a distinct physical meaning and appear to be particularly well suited 
for the description of self-gravitating fluids. 

We shall assume our system to be spherically symmetric. This is a com- 
mon assumption in the study of self-gravitating compact objects, since devi- 
ations from spherical symmetry are likely to be incidental rather than basic 
features of the process involved (see however the discussion in [27]-[29]). 

For sake of generality we shall further assume the fluid to be locally 
anisotropic (principal stresses unequal) and dissipative. 

The assumption of local anisotropy of pressure, which seems to be very 
sensible to describe matter distribution under a variety of circumstances, has 
been proved to be very useful in the study of relativistic compact objects 
and related problems (see [30] for a comprehensive review until 1997, and 
[31]-[72] and references therein, for more recent developments). 

On the other hand, it is already an established fact, that gravitational 
collapse is a highly dissipative process (see [73], [74], [75] and references 
therein). This dissipation is required to account for the very large (negative) 
binding energy of the resulting compact object (of the order of — 10 53 erg). 

In fact, it seems that the only plausible mechanism to carry away the 
bulk of the binding energy of the collapsing star, leading to a neutron star 
or black hole, is neutrino emission [76] . 
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As usual, we shall describe dissipation in two limiting cases. The first case 
(diffusion) applies whenever the mean free path of particles responsible for 
the propagation of energy is very small as compared with the typical length 
of the object. For example, for a main sequence star as the sun, the mean 
free path of photons at the center, is of the order of 2 cm. Also, the mean free 
path of trapped neutrinos in compact cores of densities about 10 12 g. cm -3 
becomes smaller than the size of the stellar core [77, 78]. Furthermore, the 
observational data collected from supernova 1987A indicates that the regime 
of radiation transport prevailing during the emission process, is closer to the 
diffusion approximation than to the free streaming limit [79] . In this case, it 
is assumed that the energy flux of radiation, as that of thermal conduction, 
is proportional to the gradient of temperature. 

The second case (free streaming), applies when the mean free path of 
particles transporting energy is larger than (or equal to) the typical length 
of the object. Since this condition may hold for a large number of stellar 
scenarios, it is advisable to include simultaneously both limiting cases of 
radiative transport, diffusion and free streaming, allowing for describing a 
wide range of situations. 

In the next section we shall introduce the relevant physical variables 
and deploy the equations for describing a dissipative self-gravitating locally 
anisotropic fluid. 

As mentioned before, a fundamental role in our study will be played by 
a set of scalars derived from the orthogonal splitting of the Riemann tensor. 
Such a splitting and the ensuing variables, are presented in Section III. In 
Section IV the physical meaning of the above mentioned scalars is discussed, 
and in Section V a general method to obtain all static anisotropic solutions 
in terms of those scalars, is presented. 

Finally the results are discussed in the last section. 

2 THE GENERAL FORMALISM 

In this section we shall present the physical variables and the relevant equa- 
tions for describing a dissipative self-gravitating locally anisotropic fluid. 
Here we shall closely follow the program outlined in [75], thus we refer the 
reader to that article for more details. 
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2.1 Einstein equations 

We consider spherically symmetric distributions of collapsing fluid, which for 

the sake of completeness we assume to be locally anisotropic, undergoing 

dissipation in the form of heat flow (diffusion limit) and/or free streaming 

radiation (free streaming limit), bounded by a spherical surface E. 

The line element is given in Schwarzschild-like (non-comoving) coordinates 

by 

ds 2 = e u dt 2 - e x dr 2 - r 2 (d9 2 + sin 2 Odcf) , (1) 

where v(t,r) and X(t,r) are functions of their arguments. We number the 

coordinates: x° = t; x 1 = r; x 2 = 9; x 3 = 0. 

The metric (1) has to satisfy Einstein field equations 

G; = kT£, (2) 
with k = Sty and which in our case read [3]: 

--n = ~ + e~\yh 0) 

-^i 1 = 4 + e " ( ^ + 7 ) ' (4) 

-kT 2 = -kT 3 = ~^(2A + A(A - v)) 

+ ^(2u" + v' 2 - AV + (5) 
4 r 

-kT w = --, (6) 

r 

where dots and primes stand for partial differentiation with respect to t and 
r, respectively In order to give physical significance to the components 
we apply the Bondi approach [3]. 

Thus, following Bondi, let us introduce purely locally Minkowski coordinates 

(r,x,y,z) 

dr = e v ! 2 dt ; dx = e x ^ 2 dr ; dy = rd9 ; dz = r sin 6d(f). 
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Then, denoting the Minkowski components of the energy tensor by a bar, we 
have 



rfiO 



?l = Tl- 



fl 



rp3 

1 3 



rp3 . 

1 3 i 



T, 



01 



r (u+X)/2 T 



01- 



Next, we suppose that when viewed by an observer moving relative to these 
coordinates with proper velocity ui in the radial direction, the physical con- 
tent of space consists of an anisotropic fluid of energy density p, radial pres- 
sure P r , tangential pressure P±, radial heat flux q and unpolarized radiation 
of energy density e traveling in the radial direction. Thus, when viewed 
by this (comoving with the fluid) observer the covariant energy-momentum 
tensor in Minkowski coordinates is 



( P + e —q — e 
— q — e P r + e 


V o o 

Then a Lorentz transformation gives 

p + P r u 2 






Pi 




\ 




p± J 



J _ J 



Tl = T\ 



2uq e(l + u) 



1 - 

Pr 



1 - ijj'^ 
2uiq 
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rp2 

1 2 



rp3 

1 3 



rp2 

1 2 



T| = -P ± 



v + \ - 



01 



e 2 T m = 



(p + P r )uje 2 



qe 2 



(1 



id 



e A 2 l 'e(l + a;) 



(7) 

(8) 
(9) 
(10) 



1 — uj 1 1 — uS 1 1 — UJ 

Note that the coordinate velocity in the (t,r,9,(fi) system, dr/dt, is related 
to uj by 



LU 



dr 
~dt 



,(A-i/)/2 



Feeding back (7-10) into (3-6), we get the field equations in the form 



p + P r 0J 2 



lujq 



6(1 



UJ 



1 

K 



+ e 



r 



(12) 
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P r + pu 2 2uoq e(l + u) 1 J 1 



+ 



+ 
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1 — UJ 2 1 — UJ 2 1 — UJ 
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1 — UT 1 — 



A + " , s • 
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Next, the four- velocity vector is defined as 



CO 



nr 



e 2 



A 

cue 2 



u 



,-,o,o), 



V(l _ w 2)l/2' (1-^2)1/2 = 

from which we can calculate the four acceleration, a a = -u^ii^ 



X-v 

uja\ = — a e 2 = 



CO 



1-cu 2 



UJUj' V* X-u ,00 X CO 



1-00' 



+ 7r) + e " 



2 I -co 2 



the shear tensor a^v and the expansion 



& liv ^/j,-!/ ~l" U fj,Q'v U u (l^ 



with 



= 9iiv ~ UfiU v 9 = u\ 



= 



-f/2 



2(l-cu 2 ) 1/2 



A + 



--X/2 



+ 



2{l-u 2 ) 1/2 



cov' + 



2co' 



+ 



0"11 



3 (1_ ^2)3/2 _ 



^(a + ^|) + 



e A / 2 U/ + 



2o/ 

1 -CJ 2 



2a; 
r 



«22 



e A r 2 (1 - a; 2 



-0"!!, 



e -V(l-a; 2 ) . 



0"33 = 



sin Ban, 



(21) 
(22) 

(23) 



(Too = ^ 2 e A eVn, 
a 0l = -cje ( ^ A)/2 o- n . 



We may write the shear tensor also as 



1 / 1 , s 



with 



(24) 
(25) 

(26) 



v , • 2o;d; s a . . 2a;' 2a; . 
e" (A + -) + e"(c;i/ + 5 ) 



1-a; 2 



1 — uj 2 r 



, (27) 



and s M being defined by 



o;e 2 



_ A 

e 2 



(l- W 2)l/2> (l_ W 2)l/2 



,0,0), 



(28) 



with the properties s M w M = 0, s M s M = —1. 

It will be convenient to write the energy-momentum tensor (7)-(10) as: 



with 



q^ = qs^, 



(29) 
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p _ P r + 2P ± 

3 ' 

p = p + e, 

Pf Pf ~\~ C ; 

Q = Q + e, 

n = p r - p ± . 

For the exterior of the fluid distribution, the spacetime is that of Vaidya, 
given by 

ds i = ^ _ du 2 + 2dudK - U 2 (d9 2 + sin 2 Odtf) , (30) 

where u is a coordinate related to the retarded time, such that u = constant 
is (asymptotically) a null cone open to the future and 1Z is a null coordinate 
= 0). 

The two coordinate systems (t, r, 9, 0) and (u, 1Z, 9, 0) are related at the 
boundary surface by 

« = t-r-Mfln(^-l), (31) 

72. = r. (32) 

In order to match smoothly the two metrics above on the boundary surface 
t = r T,{t), we require the continuity of the first and the second fundamental 
forms across that surface, yielding (see [15] for details) 

2M 

6 1 " TV (33) 
e- = 1 - ^ (34) 

[^] S = M S , (35) 
where, from now on, subscript E indicates that the quantity is evaluated on 
the boundary surface S, and (35) expresses the discontinuity of the radial 
pressure in the presence of heat flow, which is a well known result [80]. 

Eqs. (33), (34), and (35) are the necessary and sufficient conditions for a 
smooth matching of the two metrics (1) and (30) on E. 
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2.2 The Riemann and the Weyl tensor 



We know that the Riemann tensor may be expressed through the Weyl tensor 
, the Ricci tensor R a/3 and the scalar curvature R, as: 

R-aflii = Cup^ + -Rpgan ~ ^Ra/3^ + -R a ^/3 

-^^9ap-^R(S^g a ^-g a pS^). (36) 

In the spherically symmetric case, the magnetic part of the Weyl tensor 
vanishes and we can express the Weyl tensor in terms of its electric part 

(E aP = C ail35 u J u 5 ) as 



(37) 



with = g m g v p - g^g ua , and rj^ af3 denoting the Levi-Civita tensor. 

Observe that E a p may also be writen as 



with 



E a/ 3 = E(s a sp + -h a p), 



E 



A + 



A(A - v) 



AV v 1 - A' 2(1 - e A ) 



+ 



2 r 
satisfying the following properties: 

E a a = 0, E al = -E( Q7 ), EarfU 1 = 0. 



(38) 



(39) 



(40) 



2.3 The mass function and the Tolman mass 

Here we shall introduce the two most commonly used definitions for the mass 
of a sphere interior to the surface E, as well as some interesting relationships 
between them and the Weyl tensor. These will be used later to provide 
physical meaning to the five scalars quantities which will be derived from the 
orthogonal splitting of the Riemann tensor. 
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2.3.1 The mass function 

For the line element (1) the mass function m is defined by: 



p3 
-^232 



1-e 



2m 

r 



Then, using (2), (36) and (38) we may write 

3m k _ k , _. ~ x _. 



(41) 



(42) 



Observe that from (3) and (41) the mass function may also be writen as 

m = £ T r 2 T °dr. (43) 
2 Jo 



Another interesting relationship for the mass function may be obtained 
as follows. From field equations (12)-(14), and (36), (38) and (43) we get 



rn 



= ^ (T ° + Tl - T|) + 



r 3 E 



6 ^ ^ 3 

Next, differentiating (44) with respect to r and using (43) it follows 



(44) 



and integrating 



E = — 



2^3 



Finally, inserting (46) into (44) we obtain 
m(r, t) 



« r 3 T _ " 



(T )V 



6 u 6 Jo 

Now, there are three specific situations when Tq = p and T 2 2 
namely: 



(45) 

(46) 
(47) 

Tl = n, 



In the static regime, i.e. when uj as well as all time derivatives vanish. 
In the quasi-static regime, where (see [15]) 

(48) 



uj 2 ps Co ps v ps A w i>\ ps A 2 ps 
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• Immediately after the system departs from equilibrium, i.e. on a time 
scale of the order of (or smaller than) the hydrostatic time, in which 
case uj m \ m i> m 0;lu ^ 0. 

Thus in the three cases above, (46) and (47) become 

E=-£- 3 J\z(p)>dr+^n, (49) 

m(r,t) = -r 3 p - - r 3 (p)' dr. (50) 
6 6 Jo 

The first of these equations relates the Weyl tensor to two fundamental 
physical properties of the fluid distribution, namely: density inhomogeneity 
and local anisotropy of pressure. The second one expresses the mass function 
in terms of its value in the case of a homogeneous energy density distribution, 
plus the change induced by density inhomogeneity. 



2.3.2 Tolman mass 

An alternative definition to describe the energy content of a fluid sphere was 
proposed by Tolman many years ago. The Tolman mass for a spherically 
symmetric distribution of matter is given by (eq.(24) in [81]) 



m T = ^ r 2 e^' 2 (t o ° - 7? - 2T 2 2 ) dr 

+ kt^Zimr^hmV* (51) 

where L denotes the usual gravitational Lagrangian density (eq.(10) in [81]). 
Although Tolman's formula was introduced as a measure of the total energy 
of the system, with no commitment to its localization, we shall define the 
mass within a sphere of radius r, completely inside E, as 



\ Jo r2eiU+X)/2 ( T - T l - 2T l) dr 

+ I r r 2> +A )/ 2 — ( — \a a(3 dr (52) 

2 Jo dt[d[d (g°f>y/=g) /dt] J 9 d [b > 
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This extension of the global concept of energy to a local level [82] is suggested 
by the conspicuous role played by uit as the "active gravitational mass", 
which will be exhibited below. 

Now, it can be shown after some lengthy calculations (see [83] for details) 
that 



>+A)/2 



K 



m(r, t) — 2 r3 ^~i 



rriT = e 

Replacing Tj 1 by (4) and m by (41), one also finds 



(53) 



ttit 



e ("-*)/V- 



r 
~2 



(54) 



This last equation brings out the physical meaning of ttit as the active grav- 
itational mass. Indeed, as it follows from (17), the gravitational acceleration 
(a = —s v a v ) of a test particle, instantaneously at rest in a static gravitational 
field, is given by (see also [84]) 



-A/2 v , 



-v/2 



a = 



rriT 



(55) 



Another expression for which appears to be more suitable for the 
discussion in Sec. IV may be obtained as follows. Taking the r-derivative of 
(54) (see [83] for details, but notice some minor misprints in eqs.(31) and 
(38) in that reference as well as slight changes in notation) and using (44), 
(53), and field equations, we obtain 



rm' T — 3m,T 



+ 



e (v+\)/2 r 3 



e (\-u)/2 r 3 



\{Tl-Ti)-E 
A 2 Az>\ 



A + 



(56) 



which can be formally integrated to give 
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~ T l 2r ( 1 2~J (57 » 

or, using (46) 



rriT = 




For the three cases considered above (i.e. the static regime, the quasi- 
static regime, and immediately after the system departs from equilibrium) 
we have T® = p and T| — T\ = II. Thus, in any of these cases, (58) expresses 
the Tolman mass of a sphere of radius r interior to S, in terms of its value 
in the case of a homogeneous energy density and locally isotropic fluid in 
equilibrium (first term), plus the change induced by density inhomogeneity 
and local anisotropy (second term), plus changes derived from the fact that 
the system is not in equilibrium (last term). For a discussion on this last 
term see [83] .We shall, come back to this expression in Sec. IV. 

The important point to stress here is that the second integral in (57) (or 
(58)) describes the contribution of density inhomogeneity and local anisotropy 
of pressure to the Tolman mass. It is also worth noticing that when the sys- 
tem is evaluated immediately after its departure from equilibrium, the value 
of uj remains unchanged. Therefore the physical meaning of my, as the ac- 
tive gravitational mass obtained for the static (and quasi-static) case, may 
be safely extrapolated to the non-static case within that (hydrostatic) time- 
scale. 

2.4 Structure and evolution equations 

As shown in [75] the following set of equations may be derived to describe 
the self-gravitating fluid: 

p* + (p + P r )6 = 2 -{0 + |)n - $ - 2q(a + £), (59) 
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F* + (p + P r )a + — n = \{a - 46) - ?, (60) 
r 3 

6 2 a 2 + 9 2as x k . ~ ^ . . 

r + T + T ~ a ~ a = - o (P + 3Pr ) + «n, (61) 

3 6 r 2 

<2 +9) = — + ¥«■ (62) 

„ /c + 9 a* 6a as 1 a 2 . . 

B + -n = -«t_ a .____ + _ + _ (63 ) 

(f A + f> + f ) + (E - ± n 4- K -PT = (64) 
(E + 2 P " 2 n) = V ( 2 n " E) + 2 g( 2 + } ' ( } 

3/77 /■£ /■£ . ~ . ^_ 

with /t = /* = /><*, o a = as a , and 

2+* = — • ( 67 ) 

These equations are not independent and, of course, provide no more 
information than the one contained in Einstein equations, however we present 
them all, since depending on the problem under consideration, it may be more 
advantageous using one set instead of the other. 

The first two equations comes from the "conservation" equations T{? = 0. 
Equations (61) (Raychaudhuri equation) and (62) are derived from the Ricci 
identitities, whereas equation (63) is a consequence of (2), (36) and Ricci 
identities. The next two equations follow from the Bianchi identities writen 
in terms of the Weyl tensor (see [75] for details). Finally, (66) is just (42). 

We shall next present the orthogonal splitting of the Riemann tensor, and 
express it in terms of the variables considered so far. 
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3 THE ORTHOGONAL SPLITTING OF THE 
RIEMANN TENSOR 

The orthogonal splitting of the Riemann tensor was first considered by Bel 
[85], here we shall follow closely (with some changes) the notation in [86]. 
Thus following Bel, let us introduce the following tensors: 



Yap = Ra^su^u 5 , (68) 

1 

—1 



Z a p =* R ay/ 3SU J U 5 = ^VayepR^ pgrf ' U & , (69) 

Xc* =* K^su'u 5 = L oR^uiu 6 , (70) 



with R* aj3l5 - \ri epi& R a p ep '. 

It can be shown that the Riemann tensor can be expressed through these 
tensors in what is called the orthogonal splitting of the Riemann tensor (see 
[86] for details). Now, instead of using the explicit form of the splitting of 
Riemann tensor (eq.(4.6) in [86]), we shall proceed as follows. 

Using the Einstein equations we may write (36) as 

R ai ps = C a \ s + 2kT%5^\ + kT{\5^ 5] - 5 [a ^ ), (71) 
then feeding back (29) into (71) we split the Riemann tensor as 

n (38 — n (I) /3S (38 n {III) I35i 

where 

R% ps = 2k~pu^u { ^\ - 2M%5% + K (p - 3P)(|^^ - 6%6%) 
R^ p5 = 2«(n [ ^J, + qs^u [P 5\ + qu^s 



R(?n )f 38 = ^ [ Vl ] -e a \e P s,E 



(73) 



with 

e Q7/3 = u^ri^p, e aiP u p = 0, (74) 
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and where the vanishing, due to the spherical symmetry, of the magnetic 
part of the Weyl tensor (H a/ 3 =* C ai psu f u 6 ) has been used. 

From (74) it follows that e va p = u a u p r}%' yi 'rf (n , a p, producing 

e^ v e vaP = PJi% - 6£h} + u a {u% - 6^), (75) 

or, contracting a with \i in (75) 

e^e^ = -2h}. (76) 

Using the results above, we can now find the explicit expressions for the 
three tensors Y a/3 , Z af3 and X a/3 in terms of the physical variables, we obtain 

Y a(3 = jkp + 3P)h a/3 + |n a/3 + E a p, (77) 

Z a p = -qs^ea^, (78) 

and 

X a f3 = -^ph a /3 + ~j^-<*P ~ (79) 

From the above, we can obtain expressions for two quantities endowed with a 
profound physical meaning (see [86] , [87], [88] and references therein). They 
appear in the orthogonal splitting of the Bel tensor [89]. One of them is the 
Bel superenergy, defined by 

W= l - {X^X afi + Y^Y^) + Z^Z al3 , (80) 

the other is the super-Poynting vector, defined as 

P a = e aPl {Y^Z^-X]Z^). (81) 

Similar quantities may also be defined from the orthogonal splitting of 
the Bel-Robinson tensor [90]. However it should be noticed that due to 
the vanishing of the magnetic part of the Weyl tensor in the spherically 
symmetric case, the super-Poynting vector associated to the Bel-Robinson 
tensor vanishes. Also, in the conformally flat case the superenergy associated 
to the Bel-Robinson tensor vanishes. In other words, these two definitions 
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cover a much wider range of situations, when defined from the Bel tensor. 
Of course, in vacuum both set of definitions coincide. Then, from (80) and 
(81), using (77, 78, 79) we find 

W = + 6 ~ pP + 9p2 + 4n2) + l E * + Y g ~ 2 ' (82) 

K 2 

P a = ^q(P + Pr)s a . (83) 

Observe that the superenergy associated to the Bel-Robinson tensor (W), 
which is defined by 

W = E a ?E afS , (84) 

(assuming the magnetic part of the Weyl tensor vanishes, as it happens in 
our case), takes the form 

W = \E\ (85) 
Combining (85) with (82) we found 

2 2 

W -W = ^-(5p 2 + 6pP + 9P 2 + AU 2 ) + —q 2 . (86) 

It is also worth noticing that the super-Poynting vector vanishes if and only if 
there is not dissipative flux. This fact clearly illustrates the physical meaning 
of this vector, and fully justifies its name. 



3.1 Five relevant scalars 

We shall now derive five scalars quantities (hereafter referred to as structure 
scalars), in terms of which we shall write our equations (59)-(66). 

Let us first observe that tensors X af3 and Y a/3 can be splitted in terms of 
their traces and the corresponding trace-free tensor, i.e. 

X a p = -TrXh a/ 3 + X <a/3> , (87) 

with TrX = X^ and, 

X <aP> = Kh^X^ - ^TrXh^). (88) 
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From (79) we have 

TrX = X T = K p, (89) 



h a p- 



and 

X <a ^ > = X T ,is a s i f '^ r ). (90) 

where 



kY\ 

X TF = {—-E). (91) 



In a similar way we obtain 



TrY = Y T = ^(p + 3P r - 211), (92) 



h a p- 



and 

Y <a(}> = YM»n»^ (93) 

with 



Y TF = {— + E). (94) 
Finally a fifth scalar may be defined from (78) as 

Z = JZcfiZoP = -j=q. (95) 

From the above it follows that local anisotropy of pressure is determined by 
X TF and Y TF by 

kU = X TF + Y TF . (96) 

We can now rewrite equations (59)-(66) in terms of our five structure scalars 
(Xt, Xt F , Yt-i Ytf-, Z): 

*lp* + hx T +Y T +X TF +Y TF )e = he + ^){X TF +Y TF )-V2{^-+aZ+-Z), 
16 6 1 I r 

(97) 

1 1 /o r/* 

^PJ + -(X T + Y T + X TF + Y TF )a + -(X TF + Y TF ) = ^-(L-26)Z-^2—, 

Zi o V o Z Z 

(98) 

0* + l e 2 + — - at - a 2 - -as 1 = -Y T , (99) 
3 6 r 
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^ +e ^ = -2T + -r z (100) 

at + a 2 + °- + ±6><t - -/ - ^ = -Y TF , (101) 



\ [(Y T + - 2X TF + X T \ (6 + \a) + X TF y = -^yfiZ, (102) 

(|p - X TF )t = ^X TF + ^ZC- + 9), (103) 



3m 

Whereas for the Bel superenergy and the super-Poynting vector we find 



— = - X TF . (104) 



W = ±(X*+ YS) + \(X* TF + Y FF ) + Z\ (105) 

and 

a/2 

P a = -^Z(X T + Y T + X TF + Y TF )s a . (106) 

3.2 On the physical meaning of the structure scalars 

Let us now focus on the physical meaning of the scalars introduced in the 
previous subsection. 

The physical meaning of Xt and Z is evident and does not require further 
clarification. 

Let us now consider Xt F - From (103) it follows that in the absence of 
dissipation (Z — 0), (using the regular center condition), 

pt = ^X TF = 0. (107) 

In other words, in the absence of dissipation, X TF controls inhomo- 
geneities in the energy density. 

The role of density inhomogeneities in the collapse of dust [91] and in par- 
ticular in the formation of naked singularities [92]-[99], has been extensively 
discussed in the literature. 
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In the non-dissipative locally isotropic case, we obtain from (103), = 
<^ E = 0. This link between the Weyl tensor and energy density inhomo- 
geneity and the fact that tidal forces tend to make the gravitating fluid more 
inhomogeneous as the evolution proceeds, led Penrose to propose a gravi- 
tational arrow of time in terms of the Weyl tensor [100] (see also [101] and 
references therein). 

However, the fact that such a relationship is no longer valid in the pres- 
ence of local anisotropy of the pressure and/or dissipative processes, already 
discussed in [75] , explains its failure in scenarios where the above-mentioned 
factors are present [102]. 

Here we see that the single scalar Xtf (in the absence of dissipation) 
controls density inhomogeneities, and therefore should be the fundamental 
ingredient in the definition of a gravitational arrow of time. If dissipative 
processes are present, the scalar Z should be incorporated into that definition. 

To establish the physical meaning of Yt and Ytf let us get back to equa- 
tions (57), using (94) (for the three cases considered in 2.3.1) we get 



rriT 




Comparing the above expression with (58) we see that Ytf describes the 
influence of the local anisotropy of pressure and density inhomogeneity on 
the Tolman mass. It is also worth recalling that Y T f, together with X T f, 
determines the local anisotropy of the fluid distribution. 

Finally, observe that for a system in equilibrium or quasi-equilibrium, 
the Tolman mass (52) becomes 

m T = - r r 2 e^' 2 (T - 7? - 2T 2 2 ) dr, (109) 
2 Jo v ' 

which for those two regimes may be written as 

m T = r 2 e iu+x)/2 Y T dr. (110) 

Jo 

Thus Yt appears to be proportional to the Tolman mass "density" for 
systems in equilibrium or quasi-equilibrium. 
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4 ALL STATIC ANISOTROPIC SPHERES 



In this section we shall restrict ourselves to static sytems. We shall see how 
the metric corresponding to any static anisotropic sphere can be expressed in 
terms of the structure scalars. We shall explore three possible alternatives. 



4.1 First alternative 

From (41) and (104) it follows at once that 

e- x = l-^r 2 (^X T -X TF ). (Ill) 
Next, using (99) and (101) in the static case, we may write 

a=^(Y TF + Y T ). (112) 

On the other hand we have in the static case 

A V 1 . . 

a = e~-, s 1 =e~, (113) 

where (17) and (28) have been used. Feeding back (113) into (112), and 
integrating we obtain 

f 2r y tf+ y t dr 

where C is a constant of integration easily obtained from (33). 

Then, using (118) and (114), the line element (1) in the static case may 
be writen as 

f2r Y TF + Y T ^ 

ds 2 = C e 3 i-^x T -x TF) dt 2 dr 2_ r 2 ( ^2 +sin 2^ 2) 

(115) 

Thus we see that all possible spacetimes generated by static anisotropic 
fluids, are fully determined by two scalars, namely: Ytf + Yt and \X T — X TF . 
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4.2 Second alternative 

Alternatively we may proceed as follows. From (43), (89) and (104) we obtain 



r 3 / m' 



3 \ r 

which after integration produces 

X> 



m ( r ) = -\ — - X tf), (H6) 



m(r) = r 3 



(j^r + Cl ). (117) 



Or, using (41) 

1 - 2r 2 ^dr + ci) , (118) 

where the constant of integration C\ may be easily calculated from (34). 
Next, the field equation (4) in the static case may be writen as 

« ^ 1 e~ A - 1 x v' , 
or, using (41), (92), (94), (112) and (113) 

with 

Y h = l -{Y T + Y TF ). (121) 



We can now integrate (120) to obtain 



^ dr 

c 2 e , (122) 



where the constant of integration c 2 may be obtained from (33). 

Thus using (118) and (122), the line element (1) in the static case may be 

writen as 

ds 2 = cj i-^(f 4^«r +ci) dr dt 2 1 -rfr 2 -r 2 (^ 2 +sin 2 ^0 2 ), 

1 - 2r 2 (J ^dr + Cl ) V r ; ' 

(123) 

allowing the representation of all possible metrics in terms of two scalar 
functions, Y h and X TF . 
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4.3 Third alternative 



In the two previous alternatives we have seen that all line elements corre- 
sponding to an anisotropic fluid may be determined by two scalars. However 
in both cases these two scalars are a combination of four (alternative I) or 
three (alternative II) of our structure scalars. Here we shall present a third 
alternative, which allows for describing any line element in terms of only two 
structure scalars, namely Xtf and Ytf- 

Now, in the static case we may write from (39), 



E 



r v // A' 1, A' 1-e' 



then introducing new variables 



y = e 



1/ 
2~ 



u 
u 



(124) becomes 



y' + 2y 



u" - 2fl + 4 

r r* 



2u(l - 2r 2 E) 



u'- * 

r 



r 2( u /_2) ' 

which after integration yields 

y = e -J*(r)dr( J e fk(r)drf^ dr + 

with 



(124) 



(125) 



(126) 



(127) 



k(r) 



dr 



u , 



ln(u' - -) 
r 



/(0 



2u(l - 2r 2 E) 



r 2 K-f) ' 

and where Ci is a constant of integration easily obtained from the junction 
conditions. 

Then, getting back to original variables, (127) becomes 



A 



v__\ 

2 r r 



(1 - 2r 2 £) + e^r 2 Ci + e~ v r 2 J (2r 2 E)'- 



; dr. 



(128) 



Next let us introduce the new variable z by, 

e 2fzdr 



e v = 



,2 ' 



(129) 
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producing 



v 



1 1 



z(r) = - + -. 



(130) 



Using (129) and (130) in (128) we get a a link between E and z 



z\r 



) = (l-2r 2 E)+r A e-$ 2z{r)dr C l + r A e-$ 2zi - r)dr j{2r 2 E)> 

Next, from field equations (4) and (5) it follows 



3 J 2z(r)dr 



-dr. 



v" V 



2r r z 



+ e 



A' V 



2 v 2 r 



+ - - 



(131) 



(132) 



which in terms of the variables z and y introduced above, becomes 

,M 6 4 , 21 

1/' + y(— + 2z-- + —) = — - + «n . 

Integrating (133) we obtain for A: 

Z 2^ e f(^+Mr))dr 



(133) 



,A(r) 



r 6(_ 2 ; ^(r)(i+«n(^)eJ g gg rfr + c) 



(134) 



where C is a constant of integration. 

In terms of Xtf and Ytf, equations (131) and (134) may be writen as: 



(r) = l + vi [1 ~ r2{YTF ~ Xtf)] + ^ /22(r)dr (Ci + / [r 2 (Y TF - X TF )\ 



D J 2z(r)dr 



-dr, 



and 



-Mr) 



f ( ,s 4 , +2z(r))dr 

r 6 (-2 / '(Q[WTf+*W)r 2 ]e J dr + ^ 



(135) 



(136) 



Thus, given X^f and Y^f and using (136) in (135) we obtain z, which by 
virtue of (129) allows to find v. Then using z in (136) determines A. 
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This approach is essentially equivalent to the method for obtaining static 
anisotropic solutions presented in [68] (see also [103]). 

As a simple example let us consider all conformally flat anisotropic fluids. 
Conformal flatness implies Ytf — Xtf- Feeding back this condition into 
(135) we obtain for z 

o - — 

z = - + — tanh( / —dr). (137) 
r r J r 

In order to specify a single solution we have to provide another condition 
on our scalars. Thus for example, if we assume further the energy density 
to be constant, then Xtf = implying Ytf = 0, leading to the well known 
Schwarzschild interior solution. 

5 CONCLUSIONS 

We have presented a systematic study of spherically symmetric self-gravitating 
relativistic fluids, based on scalars functions derived from the orthogonal 
splitting of the Riemann tensor. In the most general case (dissipative and 
anisotropic fluid) we have five scalars, which reduce to two, in the case of dis- 
sipationless dust and static anisotropic fluids, and to one for static isotropic 
fluids. 

The motivation to present such a study and to consider further these 
scalars in the study of self-gravitating objects stems from their distinct phys- 
ical meaning. As we have seen, two of them (Xt and Z) define the energy 
density and the dissipative flux, respectively. In the absence of dissipation, 
Xtf controls the inhomogeneity of energy density and therefore is the rele- 
vant quantity in any definition of a gravitational arrow of time a la Penrose. 
Of course, in the dissipative case Z must also enter into that definition. 

The two scalars Ytf and Yr are related in a conspicuous way to the 
Tolman mass. On the one hand Ytf describes the influence of, both, energy 
density inhomogeneity and local anisotropy of pressure, on the Tolman mass. 
On the other hand Y T acts as a Tolman mass density. It is worth noticing 
that these two scalars are the only ones that appear in the "kinematical" 
equations (99) and (101). Also observe that Z is the only structure scalar 
appearing in (100). 
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In the static case, Einstein equations reduce to three ordinary differential 
equations for five variables (p, P r , P±, v, A), implying that any specific solu- 
tion is determined by two scalar functions. This was illustrated in the three 
subsections above. Furthermore, these two scalars functions may be Xtf 
and Ytf a s shown in the third alternative developped in subsection 4.3. This 
brings out further the physical relevance of the structure scalars. 
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